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ABSTRACT 

The surface abundances of extreme helium (EHc) and R Coronae Borealis (RCB) 
stars are discussed in terms of a model for their origin in the merger of a carbon- 
oxygen white dwarf with a helium white dwarf. The model is expressed as a linear 
mixture of the individual layers of both constituent white dwarfs, taking account of 
the specific evolution of each star. In developing this recipe from previous versions, 
particular attention has been given to the inter-shell abundances of the asymptotic 
giant branch star which evolved to become the carbon-oxygen white dwarf. Thus the 
surface composition of the merged star is estimated as a function of the initial mass 
and metallicity of its progenitor. The question of whether additional nucleosynthesis 
occurs during the white dwarf merger has been examined by including the results 
of recent hydrodynamical merger calculations which incorporate the major nuclear 
networks. 

The high observed abundances of carbon and oxygen must cither originate by 
dredge-up from the core of the carbon-oxygen white dwarf during a cold merger or be 
generated directly by a-burning during a hot merger. The presence of large quantities 
of 18 may be consistent with both scenarios, since a significant 18 pocket develops 
at the carbon/helium boundary in a number of our post-AGB models. 

The production of fluorine, neon and phosphorus in the AGB intcrshell propagates 
through to an overabundance at the surface of the merged stars, but generally not in 
sufficient quantity to match the observed abundances. However, the evidence for an 
AGB origin for these elements, together with near-normal abundances of magnesium, 
points to progenitor stars with initial masses in the range 1.9-3 Mq. 

There is not yet sufficient understanding of the chemical structure of CO white 
dwarfs, or of nucleosynthesis during a double white dwarf merger, to discriminate the 
origin (fossil or prompt) of all the abundance anomalies observed in EHc and RCB 
stars. Further work is required to quantify the expected yields of argon and s-process 
elements in the AGB intcrshell, and to improve the predicted yields of all elements 
from a hot merger. 

Key words: stars: chemically peculiar (helium), stars: evolution, stars: abundances, 
stars: AGB and post-AGB, white dwarfs, binaries: close 
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Extreme helium stars (EHe: spectral types O-A), RCoronae 
Borealis stars (RCB: spectral types F-G) and hydrogen- 
deficient carbon stars (HdC) are early- to late-type su- 
pergiants with atmospheres almost void o f hydrogen, but 
highly enriched in carbon (|jeffervl 120081 ; IClavtonl Il996l : 
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lAsplund et"ai]|2000l ). They display an extraordinary mix- 
ture of atomic species in ratios very different to those likely 
to have been established when the star was formed. In 
addition to hydrogen and helium, anomalies include large 
enrichments of nitrogen, oxygen, fluorine, neon and phos- 
phorus, detections of lithium and s-process elements, over- 
abundance in silicon a nd sulphur, together with a large 
range in iron abund ance (|jeffervll996l ; lRao fc Lambertll996l ; 
iPandev et ai1l2006h . It is commonly accepted that the hy- 
drogen is a remnant of the outer layers of the original star, 
the enriched nitrogen comes from CNO-processed layers, he- 
lium from CNO-processed material and partially tr iple-Q- 
proce ssed material, and the carbon from 3a reactions (| Heber] 
1 19831 ). There is increasing consensus concerning the origin of 
other elements, where the signature of material from the in- 
tershell (helium-rich) layers of asymptoti c giant branch stars 
seems unmistakable (|Pandev et alJuOOrj ). 

In order to interpret these abundances, it is necessary to 
infer some cataclysm in the history of the star by which the 
hydrogen-rich surface has been either ejected or consumed, 
whilst at the same time revealing a mixture of both CNO- 
processed helium, 3-a processed helium, and other highly- 
processed material. This same process must also have in- 
volved the creation of a cool supergi ant which is curren tly 
contracting to become a white dwarf (|Schonberner|[l98r3 ) . 

It was thought likely that EHe stars or RCB stars 
formed either following a final helium s hell flash (or late 
therma l pulse ) in a cooling white dwarf |lben et al.l 1 19831 ; 
iHerwid |2000| ), or following a merger between a carbon- 
oxygen white dwarf and a helium white dwarf in a close 
binary (|Webbinkll98i ; ISaio fc Jeffervll2002l ). Consensus now 
strongly favours the white dwarf merger origin, supported by 
evidence of evolut ion timescales, pulsa t ion masses, and sur- 
face abundances (ISaio fc Jeffer y| |2002l ; IPandev et all 120061 ; 



Table 2. Notes to Tablc[T] 



IClavt on et al. I 120071 '). Efforts are now focused on secur 
ing the abundance measurements and on the question of 
whether the merger is cold (no nucleosynthesis) or hot. 
For example, IClavton et al.l (|2007l ) argue that RCB over- 
abundances of oxygen in general and 18 O in particular are 
produced from nuc l eosynt hesis during the merger, whilst 
IPandev fc Lambertl (|2010h argue semi-quantitatively that 
no additional nucleosynthesis is required to match the ob- 
served abundances of H, He, C, N, O and Ne in EHe stars. 

Evidently, we are far from an exhaustive understanding 
of (i) the evolution and (ii) the subsequent merger of two 
white dwarfs in a close binary. In the first instance, it is 
likely that the binary will have passed through at least one 
prior phase of common-envelope evolution. In the second, 
the merger of two white dwarfs involves the total destruction 
of the less massive star and the assimilation of a subsequent 
hot disk into the survivor. Both involve non-linear processes 
on dynamical timescales. Consequently, the accuracy with 
which current surface abundances may be used to infer past 
evolution may be legitimately challenged. 

The object of this paper is to clarify the surface abun- 
dances which might arise under conservative assumptions for 
the merger of a helium white dwarf (HeWD) with a carbon- 
oxygen white dwarf (COWD) and to compare these with 
observed abundance anomalies. The question of predicted 
birthrates and galactic distribution of double-white dwarf 
mergers and their correlation with the observed distribution 
of EHe and RCB stars will be addressed in a separate paper. 
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Thus, the known situation regarding surface abundances is 
reviewed in § [2] Background assumptions and calculations 
relating to the white dwarf merger model are discussed in 
§0 particularly where these relate to the following. The mix- 
ing model used to infer merger surface abundances, together 
with the input from detailed stellar evolution and nucleosyn- 
thesis calculations, is described in §2] The inferred elemental 
abundances are discussed in § [5] Finally, the sufficiency of 
the model is discussed in terms of whether additional nucle- 
osynthesis is necessary. (§|S). 



2 THE OBSERVED SURFACE ABUNDANCES 
OF EHE AND RCB STARS 

Over the last two decades, much effort has been expended 
on the accurate measurement of the surface abundances in 
EHe and RCB stars. These have been derived primarily from 
the analysis of high-resolution spectra using model atmo- 
spheres and theoretical line profiles which assume local ther- 
modynamic, hydrostatic and radiative equilibrium (LTE). 
Model at mospheres have bee n computed using the codes 
MARCS (lAsplund et al. l ll997l ) for cool stars and STERNE 
I Jefferv et alj|200ll ) for hot stars. RCB and very cool EHe 
abundances are based on MARCS model atmospheres which 
use continuous opacities from the Opacity Project, and treat 
line opacities using opacity sampling. Published EHe abun- 
dances are based on STERNE models which use ATLAS6- 
type continuous opacities and an opacity-distribution func- 
tion com puted specifi cally for a hydrogen-poor carbon-rich 
mixture |Mollerlll990l ). 

It is important to note that although the abundances 
derived are in general consistent and not particularly sen- 
sitive to the model atmosphere structure, there are likely 
to be exceptions for individual elements. There are also 
improvements which could be made to the model atmo- 
spheres, and would lead to greater overall confidence in the 
derived abundances. For example, modern STERNE models 
include Opacity-P roject continuous opaci ties and opacity- 
sampling for lines |Behara fc Jefferv 20061). These have not 
yet been used in any detailed fine analysis of an EHe star, 
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Landscape table to go here: 
Table 1. Observed abundances for EHe and RCB stars 
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Figure 1. Observed surface abundances (log number relative to solar) versus iron abundance (same units) for extreme helium stars (blue 
squares), majority RCB stars (red diamonds), and minority RCB stars (red triangles). Upper limits are shown as arrows. The dotted line 
indicates a solar composition scaled to iron. Note that the plots for Ne, F, and Li are offset vertically; a cross indicates [Fe], [X] =0,0 
(solar abundance) in every case. 



but do demonstrate sizable differences in temperature struc- 
ture and overall flux distribution from the earlier models. 
For most RCB stars the major opacity source in the pho- 
tosphere is carbon, but the most abundant species, helium, 
is not directly observable. This has led to a "carbon prob- 
lem" whereby it i s difficult to measure t he carbon abundanc e 
unambiguously (|Asplund et al.l [l997l ; iPandev et al.1 12004 ). 



Since most EHe and RCB stars have low surface gravities, 
departures from LTE may also be important, particularly 
for the measurement of s urface gravity and other quantities 
derived from strong lines (|jefferv|| 19981 ; iPrzvbilla et aLlfeoOfj ; 
IPandev fc Lamberdlioioh . 

For the present paper, published abundance data for 
EHe and RCB stars are collated and summarised in Tables[T] 
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and [2] Abundances are cited in logarithmic unit£] such that 
€i = log n, + C, En, = 1, (1) 

where 



log T.m^i + C = log Eriio/ii + C© « 12.15, 



(2) 



and the normalisations are such that the logarithmic hydro- 
gen abundance of the Sun 



£He0 e= 12.00. 



(3) 



Measurements have been rounded to 1 decimal place 
and given without errors; the reader should refer to the pa- 
pers cited for more detail, but in general these are typi- 
cally ±0.2 — 0.3 dex. For space considerations, some ele- 
ments measured in only a few stars have been omitted. 
Likewise, the HdC stars have been omitted, since iron abun- 
dances for these are not available. The EHe stars V652 Her 
dJefferv et al.l I1999T) and HP 144941 (|Harrison fc Jeffervl 
1 19971 ; IJefferv fc Harrisonl 1 19971 ) have been excluded for a 
different reason. They are comparatively hydrogen-rich, 
nitrogen-rich and carbon-poor, suggesting a different his- 
tory. The measured iron ab undance for DYCen (5 .0) is un- 
typical of other elements (jjefferv fc Heberl 1 19931 ) ; unpub- 
lished data suggest a higher value. In this paper we assume 
an iron abundance for DY Cen scaled to that of aluminium, 
silicon and sulphur (ep c = 7.3). 

The emergence of overall patterns may be seen in Fig. \T\ 
where each panel represents a different element, different 
symbols represent different groups of stars, [X] = £j — e;© 
represents the (logarithmic) elemental abundance relative to 
the solar abundance, and [Fe] represents the iron abundance 
normalised in the same way. To understand this plot, con- 
sider that a star having the same composition as the Sun 
would appear at the origin (0,0) in every panel. Stars with 
elemental abundances scaling exactly with the iron abun- 
dance would lie on a straight line through the origin and 
having gradient unity (as indicated by a broken line). Sim- 
ilar plots have been pre sented and discussed in detail by 
Asp lund et al. | (|2000t ) and lPandev et all (|2006l ). In summary, 
their conclusions were as follows. 



2.1 Elements unaffected by evolution 

Iron: several elements appear to be representative of ini- 
tial metallicity. Fe may be adopted for spectroscopic conve- 
nience, and it is unlikely to be affe cted by H and He bu rning 
and attendant nuclear reactions. iPandev et all (|2006D find 
that Cr, Mn, and Ni vary in concert with Fe, so that these 
may also be taken as proxies for the initial metallicity. 
a— elements: Mg, Si, S, and Ca and also Ti follow the 
expecte d trend in which the abundance ratio a/Fe varies 
with Fe (|Rvde fc Lambertll2004lGoswami fc PrantzodboOOl ) 
(with the possible exception of DYCen: see above). 



Conventionally, stellar abundances are given logarithmically by 
number, normalised such the logarithmic hydrogen abundance is 
equal to 12. This convention assumes that hydrogen dominates 
the composition. In evolved mixtures, hydrogen may be vanish- 
ingly small, so the convention loses value. The formalism given 
here preserves abundance values (e^) of species unaffected, for 
example, by the conversion of hydrogen to helium. 



Aluminium: abundances follow Fe with an apparent offset of 
about 0.4 dex. 

Argon: in five out of seven EHes, Ar appears to have its 
initial abundance. 

Nickel: although Ni varies in concert with Fe in both EHes 
and RCBs, there is a disconcerting offset of about 0.5 dex 
between the two groups. At this juncture, one suspects a 
systematic error due to the use of different ions or lines in the 
two groups of stars. The difference serves as a reminder that 
caution must be exercised with all abundance measurements 
discussed here. 

Zinc: like Ni and Al, Zn correlates well with Fe, with a pos- 
itive offset of w 0.8 de?0. Again, one suspects a systematic 
error. 

Minority RCBs: [Lambert fc Raol (|l994l ) identify a subset of 
four RCBs which show lower Fe abundance and higher Si/Fe 
and S/Fe ratios than the majority. These are indicated in 
Tabled] 



2.2 Elements affected by evolution 

Hydrogen: excluding DYCen and V854Cen, the combined 
sample of EHes and RCBs have H abundances log e; in the 
range 4-8. 

Lithium: a few RCBs are notably rich in lithium, which 
must have been produced simultaneously with or subse- 
quent to the proces s which made these stars H-deficient 
(|Asplund et al.ll2000l ). 

Carbon: excluding MVSgr, the EHes show a mean carbon 
abundance log = 9.3, and a range from 8.9 - 9.7, corre- 
sponding to a mean C/He ratio of 0.006 and a range from 
0.003 to 0.010. The carbon abundance is more difficult to 
measure reliably in RCBs; the mean indicated by Table [1] 
is apparently lower than in the EHes; this is probably a di- 
rect consequence of the carbon problem referred to above 
(|Asplund et al.lll997l ). 

In RCB and HdC stars cool enough to show CO, the 
12 C/ 13 C ratio is generally found to be greater than lOfJf] 
(|Warnerlll967l ; ICottrell fc Lambert! fl982l ). confirming a 3-a 
or helium-burning origin for the carbon excess. 
Nitrogen: nitrogen is enriched in the great majority of EHes 
and RCBs abo ve that expec ted according to the Fe abun- 
dance (Fig. [[]). iHeberi (|l983l ) and subsequent authors point 
out that the N abundances in general follow the trend ex- 
pected by the almost complete conversion of the initial C, 
N, and O to N via the H-burning CN and ON cycles. The 
exceptions are again DY Cen (very N-rich for its Fe abun- 
dance), and LSS99 (very little N enrichment). 
Oxygen: abundances relative to Fe range from underabun- 
dant by more than 1 dex to overabundant by almost 2 dex. 
The stars fall into two groups. Six EHes and a comparable 
number of RCBs with [0/Fe]> 1 stand apart from the re- 
mainder which have an O abundance closer to their initial 
value. The O/N ratio for this remainder is approximately 
constant at O/N « 1 and independent of Fe. 

Both groups present problems. There is no obvious 



2 At very low metallicity, Zn is thought to be enhanced by the 
s-process, but not at levels which concern us her^ 

3 the exception being VCrA with 12 C/ 13 C ~ 4 iRao fc Lambert! 
|2008|) 
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means to produce [0/Fe]> 1 (in most of these cases it is Table 3 Galactic mcrgcr ratcs for doublc whitc dwarf binaries, 

impossible to distinguish 16 O from 18 O). For the remainder, 



nearly all have an N abundance indicating total conversion 
of initial C, N, and O to N via the CNO cycles, so that an 
observed O abundance so close to the initial abundance is 
unexpected. It could be partially accounted for by dredge- 
up of 16 O from the CO-core during post-AGB evolution, as 
has b een suggested to explain high O abundances in PG1159 
stars i| Werner fe HerwidboOrj ). 

Another solution is suggested by the remarkable discov- 
ery that the 18 0/ 16 ratio in RCB stars (where observed) 
is close to and sometimes greater than unity, a ratio many - 
hund reds of times higher than expected (|Clavton et al.ll2005l . 
120071). In HdC stars, the 18 Q/ 16 Q ratios are even higher 



1 Garcfa-Hernandez et all 12009] . l20ich . However, there is as 
yet insufficient evidence to indicate that the excess O is in 
the form of ls O in all RCBs or in any EHes. 

Knowing that early-asymptotic giant branch stars pos- 
sess a high 18 Q/ 16 Q pocket in a thin layer of the He-burning 
shell, IWarnerl (|l967h speculated that under unusual circum- 
stances a star might strip all of its e nvelope material pre - 
cisely down to this narrow layer, but ICIavton et al.l (|2005l ) 
concludes this to be highly improbable. 
Fluorine: the discovery of very substantial quantities of flu- 
orine, first in several EHes and subsequently in most RC Bs, 
was also unexpected (|Pandevll200r3 ; IPandev et alj|2008h . It 
appears to be uncorrelated with Fe or O and is overabun- 
dant by 2 - 4 dex. F is produced in the He-intershell of an 
AGB star through a (complex) combina tion of a—, n- and 
p-capture reactions |Lugaro et al.ll2004l ). a conclusion con- 
firmed by an observed corr elation between C and F in the 
atmo spheres of AGB stars (|jorissen et al.|[l992l ; lAbia et al.l 
|2010| ). and by observations of post-AGB star s that show F 
at th e level predicted to be in the intershell (jWerner et al.l 
I2009D . 

Neon: high overabundances derived from Nel lines for a 
few intermediate temperature EHes were originally treated 
with scepticism - non-LTE being a possible culprit. A sim- 
ilar overabun dance measured from Nell lines in LS IV+6°2 
|jeffervl 1 19981) effectively substantiate d the Nel results in 
other stars. Pandev fe Lambertl (|2010l ) have made a recent 
non-LTE analysis of the neon abundances in three EHes, and 
confirmed a substantial overabundance approximately inde- 
pendent of the star's iron abundance. 22 Ne is produced via 
two a-captures on 14 N, so should be abundant in carbon-rich 
material derived from helium produced by the CNO-process. 
Phosphorus: an ove rabundance of P was fi r st re marked 
in BD+10°2179 by iHunger fe Klinglesmithl l| 19691 ) This 
was d iscoun ted from ultr a violet spectroscopy by iHeberl 
{ 19831 ) and iPandev et alj (|2006h ■ Overabundances have 
been reported in sever al oth er EHes by int e r alia 
Kaufmann fe Schonbernerl (|l977l ); Ijefferv fe Heberl j 19921 . 
19931 ); Ijefiervl (|l993l ). where they are systematic ally larger 
than in the sample studied bv lPandev et al.l (|2006l ). Whether 
this represents a problem with p /-values deserves further in- 
vestigation. P can be produced through neutron captures in 
an asymptotic giant branch star (of which more later). As 
the observations stand, P overabundances may be a key di- 
agnostic of previous history. 

Heavy elements: two EHes are severely enriched in Y and Zr: 
V1920 Cy g and LSE 78 with o verabundances of about a fac- 
tor of 50 jPandev et alll2004l ). A third, PVTel, is enriched 



Source 


He+He 


He+CO : CO+He 
yr- 1 


CO+CO 


Webbink (1984) 1 


2.9 


1.9 


1.2 






Xl0 _11 pc- 2 yr- 1 




Iben et al. (1996) 




0.0023 




Han (1998) 2 


0.0112 


0.0154 


0.0044 


Nelemans et al. (2001) 




0.0044 




Yu & Jefferv C2010) 




0.0027 





It is not clear whether Webbink reports the double-degenerate 
birth rate or merger rate. 

2 Model Set 1: merger rate = birth rate assumes all DD's will 
merge 

by a factor of about 10. Five other stars for which mea- 
surements were possible are considered to have their initial 
abundances of Y and Zr. Y and Zr overabundances are at- 
tributed to contamination by s-process products. The origin 
of these has not been identified. 

Only upper limits are reported for rare earth elements 
La, Ce, and Nd, all consistent with the observed abundances 
of Zr and Y. The cool EHe LSIV-14°109 has a Ba abun- 
dance consistent with its initial metallicity. 



2.3 Key Questions 

The surfaces of RCBs and EHes primarily exhibit CNO- 
processed helium. In addition, they show contamination by 
a residue of hydrogen, by 3a-processed carbon, and by ad- 
ditional ^-capture products. 

The primary challenge is to demonstrate a mechanism 
which will deliver the observed mixture, or range of mix- 
tures, at the stellar surface. This mechanism must also be 
able to explain large overabundances of Li, ls O, 19 F, possi- 
bly P, and various s-process elements. 



3 MODELS FOR THE EVOLUTION OF 
CO+HE WD MERGERS 

While evidence has accrued in favour of an origin involv- 
ing the merger of a carbon-oxygen white dwarf with a he- 
lium white dwarf, this has not always been the favoured 
model. Valid questions include whether such mergers can 
occur, with what frequency, and with what outcomes. 



3.1 The formation of CO+He binary white dwarfs 

IWebbinkl l|l984h first recognised that one consequence of 
close binary evolution would be the formation and evolution 
of double white-dwarf binary (DWD) systems that could ul- 
timately merge. 

The evolution of a main-sequence star results in an ex- 
pansion that will bring it into contact with a sufficiently 
nearby companion. Such contact may result in stable Roche 
lobe overflow, or dynamical mass transfer, resulting in the 
formation and ejection of a common envelope. The out- 
comes depend on the binary mass ratio a nd on the struc- 
ture of the larger star, and are diverse (|lben fe Tutukovl 
Il985l) . The remnants include double helium white dwarfs; 
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q=5/6 +f(m) a /2 . 



q=2/3 



q = 5/6 +«m) LB /2 



Figure 2. Stability limits for mass transfer in close double 
white dwarfs. Above the appropriate limit, mass transfer should 
be dynamically unstable. The li mits $lb ar *d $ZS refer to cold 
WD mass-radiu s relati ons due to lLvnden-Bell fc Tou tl feOOll) and 
iNelemans et al.l i200lh respectively. 



carbon-oxygen plus h elium white dwa rfs, and double carbon- 
oxygen white dwarfs. IWebbinkl (f 19841 ) estimated birth rates 
for the formation of each of these systems; similar estimates 
have formed one output of many subsequent binary-star 
population-synthesis studies of the Galaxy. 

Up to the mid 1980's, a criticism of theory was the 
absence of hard observational evidence that short-period 
white-dwarf binaries actually do form, whether as a con- 
sequence of close-binary evolution or otherwise. This prob- 
lem was largely addressed by the discovery of sign i ficant 
numbers of such systems (|Saffer et al.l 1 198 j; iMarshl [l995 
iMarsh et al.ll 19951 ) . Further discover ies were made as a result 
of large-scale white dwarf surveys dNapiwotzki et all I20031 : 
INelemans et al.ll2005l ; iMorales-Rueda et al.ll2005h . 

At present, there exists a qualitative agreement between 
observed DWD space densities and their predicted birth 
rates. The question, as it applies to DWD mergers, will be 
addressed in more detail elsewhere (Jeffery et al. in prepa- 
ration) . 



The second principle is that if the mass ratio 
q = m 2 /mi > g cri t = g H ^ — ' 



(5) 



the increase of radius due to the reduction of mass (£(m) = 
dlnr/dln?n) will exceed the increase in the Roche ra- 
dius caused by the transfer of angular momentum. Mass 
transfer then becomes dynamically u nstable and proba- 
bly causes the components to coa lesce (jPringle fc Webbinkl 
1 19751 ; iTutukov fe Yungelsonlll979r ). For this paper we have 
adopted t he mass-radius relation for a cold white dwarf re- 
ported bv lLvnden-Bell fc Tout] (|200ll ), with /3 = 1.137 and 
He = 2.02. Note that this relation, which is valid from sub- 
planetary masses through to relativistic wh ite dwarfs, gives 
£(m) markedly different to that used by INelemans et al.l 
(|200lf ). which lies very close to the classical non-relativistic 
r oc mr 1//3 relation (Fig. [2}. 

For smaller q, mass transfer will be stable; but if the 
mass-transfer rate exceeds the Eddington rate, the enve- 
lope of the accretor will heat and expand leading to the 
possible formation of a comm on envelope and which could 
also cause the stars to coalesce l|Han fc Webbinkl 1999T I . Only 
at the most extreme mass ratios will mass transfer be sta- 
ble, possibly leadi ng to the formation of AM CVn systems 
(He+He WDBs) (|Nather et al.lll98ll ). Significantly, all of 
the DWD systems fo r which mass ratios could be measured 
jMaxted et al.ll2002l ) have q > q CI it- 

Using this information together with DWD birth rates, 
population synthesis studies indicate a DWD merge r fre- 
quency for the Galaxy of between 2.3 — 4.410 ~ 3 y~ 1 jibed 
1 19901 ; INelemans et al. f boOll ; lYu fc Jeffervl I2010T ) . Other esti- 
mates are indicated in Table [3] broken down by binary type 
wherever possible. 

Of particular interest, of course, is the frequen cy of dou 



1984 




2004 





Iben fc Tutukovll 19841 ) but arguable (|Saio fc Nomotd 
Yoon fc LangerN2005h source of Ty pe la supernovae 



(SN l a). Other outcomes a r e mor e likely (|lben fc Tutukovl 
1 19841 ; ISaio fc Nomptolll995l,Jl99Sl) . Doubl e HeWD m ergers 



may produce sdO jWebbinklfl984l ') or sdB (|lbeni ri990) stars, 
which are ubiquitous in old stellar populations ([Brown et al.l 
l200ll ; iBusso et al.ll2005l ; iLee et al.|[2005l ; iRich et al.ll2005h . 



3.2 Gravitational- wave radiation and dynamical 
mergers 

There are two principles behind the idea that close-binary 
white dwarfs will merge to form a single star. 

The first principle is that angular momentum is re- 
moved from the binary by means of gravitational radiation 
(GR) and that, within a Hubble timescale, the less massive 
and consequently larger white dwarf will eventually fill its 
Roche lobe and a phase of mass transfer will begin. The 
timescale for orbital decay by GR is given by 



(r/y) = lO 7 (P/h) 8/3 M _1 (M/M )^ 2/3 , 



(4) 



where P is the orbital period, M — mi+ Bi2 is the total mass 
of the system and u is t he reduced mass ijLandau fc Lifshit3 
ll958l : lMarsh et al.lll995h . indicating that DWD systems with 
P < a few hours will reach contact within a Hubble time. 



3.3 SPH simulations of the dynamical merger 

Several simulations of the white dwarf merger process 
hav e been attempted using smoothed particle hydrodynam- 
ics (Bcnz et al. 1990: ScErctain et al. 1997: Guerrero ct al. 
120041 ; lYoon et all 120071 ; lLoren-Aguilar et all 120091 '). For 
CO+He WD mergers (e.g. 0.6+0.4 M Q ), these demonstrate 
the total disruption of the low-mass WD within roughly one 
orbital revolution (~ 90s) and the conservation of ~ 99% 
of its mass within a thick Keplerian disk. They also demon- 
strate substantial prompt heating of the disrupted material, 
with temperatures m omentarily reaching several 10 9 K 111 
the equatorial plane (|Guerrero et al.l |2004|). However, the 
temperatures are not extremely high, the degeneracy of the 
disrupted material is rapidly lifted, and any thermonuclear 
activity is ultimately quenched. 
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3.4 Nucleosynthesis during a dynamical merger 

In previous discus sion of the post-merger product, 
ISaio fc Jeffervl l|2002l ) made the simplifying assumption that 
no nucleosynthesis occurs during the merger - this is the 
cold merger approximation. 

Where 12 C and 4 He mix at sufficiently high temper- 
atures, some thermonuclear activity will occur. Any 14 N 
will also be briefly exposed to a-burning. Using an elegant 
one-zone model i n which orbita l energ y is converted to heat 
in a debris disk, IClavton et all <|2007h showed that certain 
nuclear abundances could be demonstrably altered during 
the merger. In particular, surplus ls O could be produced 
through prompt nucleosynthesis of 14 N and 4 He from the de- 
bris of the helium white dwarQ, without being subsequently 
dest royed by an additional q-ca pture to form 22 Ne. 

lLoren-Aguilar et al.l |2009l ) included a limited nuclear 
network in their dynamical merger simulation, and reported 
nuclear yields for various DWD progenitor combinations. 
Models which include some nucleosynthesis during white 
dwarf destruction will be referred to as the hot merger ap- 
proximation. 

Since the physics of DWD mergers is of substantial 
wider interest for the production of hot subdwarfs and SN 
la, the question of whether mergers are hot or cold is par- 
ticularly relevant. Rephrasing: to what extent does nucle- 
osynthesis occur as a direct consequence of heat generated 
by orbital energy dissipated during the merger, and do any 
nuclear products play a role in the subsequent evolution? 
In particular, is a hot merger necessary to explain the high 
ls O/ 16 ratio observed in RCB stars? 



3.5 Models of thermal and nuclear evolution after 
a merger 

The evolution of a WD rapidly accreting helium was first 
considered long before the possibility of DWD m ergers 
was widely recognised (|Nomoto fc Sugimotol Il977l ). This 
and subsequent calculations pursued the evolution of the 
accretor through and beyond off-centre helium ignition 



( Nomoto 



Hashim oto 198 7|; iKawai et al.|[l987l 
Jcffcrv 2000, 



ll99d ; ISaio fc No moto 199ci ISaio 



1988i; llbenl 



20021 ). 



Such models have been used to approximate evolution 
following a dynamical merger by making some working as- 
sumptions. These include the less massive white dwarf be- 
ing completely disrupted by the merger, forming a Keplerian 
disk and subsequently being assimilated onto the surface of 
the accretor. Assimilation has been assumed to be by spher- 
ical accretion at half the Eddington rat(jf| (~ 10 -5 Moy -1 ). 
Accretion was switched off once a pre-selected final mass 
was attained. 

Conceptually, this approach is flawed. It implies that, 
following shell-helium ignition and subsequent expansion, 
the reservoir of material to be accreted (i.e., the rem- 
nant of the disrupted white dwarf) remains in a Keplerian 
disk deeply embedded within the giant envelope. Although 



4 The reactions involved are: 14 N(o, 7) 18 F(/3+) ls O 

5 This rate was chosen to avoid a runaway explosion (for low ac- 
cretion rates) and to satisfy energy conservation. Higher accretion 
rates would be possible if heat could be removed asphcrically. 



such a disk might survive, it runs counter to the princi- 
le that viscous disks collapse on a much shorter timescale 
Lvnden-Bell fc Pringlel Il974l ). The disk is more likely to 
heat and expand t o form a high-ent ropy envelope in hydro- 
static equilibrium (|Yoon et alj|2007h . 

In their simulations ISaio fc Jeffervl l|2002ft assumed all 
accreted material to be deposited on the surface of the accre- 
tor, and to have a composition defined by the mean compo- 
sition of a helium white dwarf. The chemical structure of the 
accreting CO white dwarf was obtained by evolving a star 
from the zero-age main sequence through to an appropriate 
point on the white-dwarf cooling sequence. 

During post-merger evolution, there are two phases of 
convective mixing. The first is bottom-up nuclear-driven 
con vection immedi a tely followi ng shell-helium ignition. In 
the ISaio fc Jeffervl (|2000l . |2002| ) models, this occurs when 
the envelope mass is small, and does not produce very much 
mixing. Subsequent top-down opacity-driven convection de- 
velops when the star becomes a giant, and does not reach 
layers enriched by post-merger nucleosynthesis. The mod- 
els consequently show very little chemical enrichment at the 
surface from C, O or other nuclear products due to this mix- 
ing. 

The relative absence of 12 C or 16 O in thel aio fc Jeffervl 



|2002l ) COWD intershell contradicts the substantial en- 
hancements predicted in the AGB intershell (e.g. iHerwid 
(|2000l )). possibly because chemical evolution through the 
thermal-pulsing AGB was not treated in sufficient detail. 
Some additional merger sequences were therefore computed 
with enhanced /3i2 C = 0.2 and /3i6 = 0.05, where /3, repre- 
sents the mass fraction of species i. 

Consequentlv. lSaio fc Jeffervl (|2002h argued that the dy- 
namical merger would have to disrupt the outer layers of the 
COWD in order to explain the observed EHe and RCB sur- 
face abund ances of C. Scrutiny s hows that th e simple recipe 
adopted bvlSaio fc Jeffervl (|2002|) , and a lso bv lPandev et al.l 



Pandev fc Lambertl (|2010h . requires further re- 



|200ll . l2006 

finement. 

First, the COWD models would benefit from a more 
realistic abundance distribution, particularly in the helium 
layer, which should have the composition of the intershell 
region of the progenitor AGB star. 

Sec ond, it was not apprec iated that some COWD mod- 
els (e.g. ISaio fc Jeffervl I2002T ) might contain a substantial 
pocket of 18 O at the interface between the CO-core and the 
He intershell, as well as a reservoir of 22 Ne in the CO-core. 
Thus, if the outer layers of the CO-core are disrupted during 
the merger, substantial 18 and 22 Ne will be dredged up in 
addition to 12 C. 

This question is somewhat open. Not all post-AGB 
models show this 18 O pocket. Depending on the inter- 
shell temperatures, 14 N may be completely destroyed by a- 
capture to 22 Ne before the helium-burning shell passes it into 
the CO core. 

The object of the following sections is therefore to refine 
the simple recipe used in previous discussions of EHe and 
RCB surface abundances, to incorporate a more realistic 
description of the chemistries involved, and hence to develop 
a more quantitative framework in which to discuss white 
dwarf mergers as possible progenitors. 
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Figure 3. Masses of the merger components as a function of the initial mass of the COWD progenitor (m agb ) and of the initial system 
metallicity [Z = 0.0001,0.004,0.008 and 0.02). rriHc:CO represents the final mass of the AGB intershell before the star contracts to 
become a COWD. mcO:CO represents the mass of the carbon-oxygen core at the same time, mjie represents the mass of the HeWD 
secondary given by Eq. 1121 rri{ represents the total mass of the product, assuming a conservative merger. 



4 THE SIMPLE RECIPE 



4.1 Masses 



As is evident from the preceding summary, the processes by 
which two white dwarfs may merge are far from straight- 
forward. In addition, the overall distribution of elements 
in both stars is a function of the initial masses, metallic- 
ity and period in the original binary system. To become a 
COWD, the ultimately more massive star must first ascend 
the asymptotic giant branch, during which time the helium- 
rich intershell will be processed by a series of thermal pulses, 
whilst the hydrogen-rich envelope will be enriched by multi- 
ple convective dredge-up episodes. To create a close binary 
which will ultimately merge, the system must pass through 
one or more common-envelope phases, producing an abrupt 
change in the mass (at least) of one or both components. 
During the merger, the less massive white dwarf will be 
completely di srupted to form (a) a Keplerian disk and (b ) 
a hot corona jYoon et al.l 120071 ; lLoren-Aguilar et all [2003 ). 
During the disruption episode, temperatures and densities 
may become high enough for nucleosynthesis. Shear mixing 
between the disk, corona and the surface of the more mas- 
sive white dwarf seems inevitable, though how deep such 
mixing would be remains to be explored. Subsequently, ma- 
terial from the disk and corona will be accreted onto the 
surface of the more massive white dwarf; the assumption 
is that it will be chemically homogenised. As this hot ma- 
terial forces the star to expand and cool, surface convec- 
tion zones will develop from the surface. Flash-driven con- 
vection will develop following helium-shell ignition. Most 
of these processes are poorly understood. Where numerical 



models do exist, e.g. for the dy namical phases (|Yoon et al.l 
120071; lLoren-Aguilar et al.ll2009l ), or for the nuclear phases 
( baio fc Jeffervlhooi r they are not yet joined up. 

Until such time as they are, some simplifying assertions 
allow us to make order-of-magnitude arguments. The prin- 
cipal of these is that all material from the helium white 
dwarf, and all hydrogen-rich and helium-rich material from 
the carbon-oxygen white dwarf will be fully mixed during 
the dynamical phase of the merger. From the point of view 
of developing a recipe for calculating the chemical signature 
of the merged product, this assertion gives us the first set of 
parameters that will be required, namely the masses associ- 
ated with each layer of material to be mixed. We adopt the 
not ation m, : fc to represe nt the mass of layer j of star k (see 
also lSaio fc Jeffervlf2002l '), thus: 

JTiH:He: the mass of the hydrogen-rich surface layer of the 
HeWD; 

fiHe:H e : the mass of the helium core of the HeWD; 
mH:Co: the mass of the hydrogen-rich surface layer of the 
COWD; and 

wiHe:Co: the mass of the helium shell of the COWD. 

This principal assertion actually takes two forms. In the con- 
servative case, all material from both white dwarfs is in- 
cluded in the merged product. In the non-conservative case, 
some mass may be lost under which circumstance the above 
parameters represent the layer masses which survive the 
merger. SPH calc ulations suggest that white-dwarf mergers 
are c onservative (| Guerrero et al.l |2004| ; lLoren-Aguilar et al.l 
120091 ). 

Without additional processing, none of these layers con- 



10 C.S.Jeffery, A.I.Karakas & H.Saio 



tains sufficient carbon to account for the abundances ob- 
served in EHe, RCB and HdC stars. As in previous applica- 
tions of this recipe, it is necessary to make a second assertion 
that some material from the outer edge of the carbon-oxygen 
core of the COWD has been somehow included into the mix- 
ture. For now, we define: 

mco-.co'- the mass of the carbon-oxygen core of the 
COWD; and 

In addition, we have the total masses: 

mH C : the total mass of the HeWD; 

mco- the total mass of the COWD; 

mf. the final mass of the merged product. 

In the conservative case, the number of variables is re- 
duced since: 



mf 



m He + mco, 



(6) 



m k = Zjm j:k ; j = H, He, CO, k = He, CO. (7) 

Also, since helium (in general) dominates the final mixture 
and 



mHe:Hc >> "lHc:CO >> mH:CO ~ m,H:He, 



(8) 



the observed hydrogen abundance naturally constrains the 
masses of the hydrogen-rich layers: 



Ph/Phc « (m,H:CO + I«H:Ee)/mHe:He 



(9) 



where pi represents the abundance by mass fraction of 
species i (Ej/3j = 1). Since in general EHes and RCBs show 
eh/ehc « Pn/ipKe/A) < 10~ 4 , we conclude 



m,H:CO + m,H:Hc < 10 mHc:Hc 



(10) 



We have attempted to use a similar argument to esti- 
mate the mass of the carbon-oxygen core to be included in 
the mixed layers (m m i x: co) by supposing that the observed 
C and O come entirely from the COWD core; 



(/?C + /9o)//3He ~ m m i X :Co/W-He:He 



(11) 



It will be seen that this approach is too crude. 

With five masses to be adjusted in order to account for 
the abundances of hydrogen, helium, carbon, nitrogen and 
oxygen, the recipe appears under-constrained. Fortunately, 
stellar evolution t heory provides a dditio nal information. 

For example, ISaio fc Jeffervl (|2002h noted that in the 
conservative case for a 0.6 Mq carbon-oxygen white dwarf 
merging to form a 0.9 Mg product, the donor would likely 
have been predominantly helium with a mass mHc:Hc = 
0.3 Mq. On its surface would have been a hydrogen-rich 
envel ope of mass mne:H ~ 10" 3 - 1O" 4 M (|Driebe et all 
1 19981 ). reduced to this value by Roche lobe overflow during 
its first ascent of the giant branch. 

We now introduce the notion that a similar connection 
exists between the carbon-oxygen white dwarf and the initial 
star in the binary system. Assuming that both stars evolve 
as single stars up to the point of merged, then models of 

6 The fully self- consistent approach would be to find the initial 
binary (mi,rri2,P rb) that will produce a close white dwarf pair 
of appropriate dimensions and then to compute the evolution of 
both components in detail, including their passage through any 
mass transfer or common-envelope phases. 



stellar evolutio n through to the late asymptotic giant branch 
(|Karakasll20iOh (Fig. [3j) give values for mHe:CO and mco:CO 
in terms of 

m a gt>: the initial mass of the star which becomes a COWD. 

Such models also give a value for the mass of the AGB star 
hydrogen envelope, but most of this will be substantially 
removed by stellar winds and may form a planetary nebula 
before the star becomes a white dwarf. 

Establishing the mass of the COWD as a function of 
its progenitor mass, the mass of the HeWD must lie below 
the minimum mass for core helium ignition, approximately 
0.48 Mq , and above the critical value for stable mass transfer 
(§ I3,2|l , To restrict the number of free parameters in our 
model, we therefore set 



mHe = qtritmco, 



(12) 



the lower limit for dynamical mergers (Eq. [5J- This auto- 
matically prescribes the ratio mHc:Hc : mH c: co for a given 
mco and Z, and hence determines the dilution by the he- 
lium white dwarf of elements produced in the AGB inter- 
shell. However, we note that more massive HeWDs may ex- 
ist, and that less massive HeWDs may merge as a result of 
a common-envelope phase. 

Thus, our simple recipe for predicting the surface abun- 
dances of the product of a cold, conservative He+COWD 
merger now only requires mi : co and Z as primary inputs, 
together with a prescription for the composition of each com- 
ponent of the mixture. 



4.2 Composition 

We introduce the notation j3ijk to refer to the mass fraction 
of species i in layer j of component k. Generally i = z, the 
atomic number. We currently consider the abundances of: 
4 H, 4 He, 12 C, 13 C, 14 N, 16 0, 18 0. 19 F, 22 Ne, 23 Na, 24 Mg, 
27 Al, 28 Si, 31 P, 32 S, 40 Ar, 40 Ca, 48 Ti, 51 V, 52 Cr, 55 Mn, 56 Fe, 
59 Co, and 59 Ni. 

Assuming the binary system was established with an 
initial metallicity Z, the composition of the hydrogen-rich 
layers in both components is then defined by a scaled solar 
composition (/Si©). We adopt: 



/3i:H:He = <*i {Z)f3iQ . Z/Zq , 
/^HoiHiHc = 0.28, 
/3h:H:Hc = 1 — £i> 1 /3i : H:He, 
/3i:H:CO = /3i:H:H C . 



i > 2, 



(13) 
(14) 
(15) 
(16) 



Generally, a; = 1, but in metal-poor environments 
(Z < Z Q /10) the abundances of ie O, 18 0, Mg, Si, S, 
Ca, Ti and Mn are observed to exceed the scaled solar 
value by as much as 0.5 dex (|Goswami fe Prantzosl l2000l : 
iRvde fc Lambertll2004r ). cti(Z) has been chosen accordingly. 
The last relation implies that dredge-up into the H-envelope 
has been ignored, especially dredge-up on the AGB, but also 
dredge-up prior to the AGB. This is justified because mH : co 
is so small that, apart from hydrogen, this layer makes a 
negligible contribution to the merger composition. 
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4.2.1 He WD core 



4.2.3 COWD core boundary 



The composition of the He WD core is assumed to have been 
produced by CNO-cycle hydrogen burning. This converts 
practically all of the carbon and oxygen (depending on tem- 
perature) to 14 7V. To allow for incomplete CNO cycling in 
low-mass stars, we introduce the branching ratio /cno be- 
tween the full CNO cycle and the CN cycle. The composition 
is then given by: 

ftiHeiHc = /3i : H:Ho, except . . . (17) 

^12C:Hc:Ho = 0, (18) 

/3l3 C:H c:Ho = 0, (19) 

/3l4 N;Hc;Hc = Ec,N/3i:H:Ho + /CNO £o/3i:H:He, (20) 

^i6 0:Ho;Hc = (1 - /cN0)/3l6 :H:He, (21) 

^OiHetHe = (1 - /CN0)/Jl8 :H:Hc, (22) 

/3H:Hc:Hc = 0, (23) 

/3He:Hc:Hc = 1 — £i^2,3i:He:Hc ■ (24) 

Ei represents a sum over all isotopes of species i. So far, we 
have always used /cno = 1, implying the full CNO cycle. 



4.2.2 COWD shell 

The helium-rich layer of the COWD corresponds to the in- 
tershell of the progenitor AGB star. The composition of this 
layer is the most interesting of all since it contains a combi- 
nation of CNO-cycled helium, various a-capture products, 
and the products of a nuclear network which includes s- 
process neutron-capture products. The yield of each isotope 
from these processes is a sensitive function of the initial 
abundances and of the temperatures and densities through- 
out successive thermal-pulse cycles of the AGB star. Given 
the importance of this layer to the final composition of the 
merged product, we have adopted intershell com positions 
from a grid of full AGB-star evolution calculations (|Karakasl 
|2010| ). These provide fractional abundances for 4 He, 12 C, 
le O, 17 0, ls O, 19 F, 20 Ne, 21 Ne, 22 Ne, 23 Na, 24 Mg, 25 Mg, 
26 Mg, 27 Al, 28 Si, 29 Si, 30 Si, 31 P, 32 S, 33 S, and 34 S and define 
/3i:He:CO for these species. In our recipe, the isotopes of Ne, 
Mg, Si and S, are combined, since these are not resolved 
observationally. 13 C and 14 N are destroyed in the intershell, 



with 14 N being converted to 22 Ne. Thus: 



/?H:Hc:CO = 0, 
/3l 3 C:Hc:CO — 0, 
/3 14 N:He:CO ~ 0. 



(25) 
(26) 
(27) 



The AGB model grid of [Karakaj (j2010T ) provides informa- 
tion about light-element abundances, core and shell masses 
in post-AGB stars over a large range of m ag b and Z, but 
lacks detailed information for s-process elements. Consistent 
calculations for Zn, Y, Zr, Ba, La, Ce and Nd are available 
for one model with m ag b = 2 Mq, Z — 0.0001 (Lugaro et al. 
2011, in preparation) and one with m ag b = 3 Mg, Z = 0.02. 
We have used these as indicative of the range of s-process 
yields in AGB stars of intermediate initial mass. 



The outer layers of the carbon-oxygen core of the COWD 
will obviously lack hydrogen and helium and be dominated 
by 12 C and 16 from the 3-a and 12 C(a,7) 16 reactions. 

An initial approach was to adopt a factor fco = 
/3i2 C /(/3i2 C + /3i6q) = 0.8 based on typical abundances for 
COWD cores from contemporary computations. Since the 
outer layers of the core will be those which most recently 
exited the base of the intershell, the abundances of most 
other elements can be set equal to those of the intershell: 



^H:CO:CO = 0, (28) 
/3 He :CO:CO = 0, (29) 

Pl6 . C0 . C0 = (1 - /co)(^H 0: Hc:CO + /3l2 C :Ho : CO + ^180:He:Co) 

(30) 

/5l2C:CO:CO = /co(/5He:He:CO + /3l2C:Ho:CO + /3l60:Hc:Co) 

(31) 

ft:CO:CO = A:Hc:CO, %± 1,2,6,8 (32) 

Apart from the introduction of realistic intershell 
abundances, this a pproach follows that discussed by 
ISaio fc Jeffervl (|2002l ), which failed to explain the oxygen 
abundances without recourse to an ad hoc argument. Follow- 
ing dis coveries of a large 18 O excess in several RCB and HdC 
stars (jClavton et al.l 120071 : iGarci'a-Hernandez et al.| 20091 ) . 
the models for COWDs used bv lSaio fc Jeffervl (|2002h were 
re-examined and found to contain a substantial pocket of 
18 O at the interface between the CO-core and the He inter- 
shell, as well as a reservoir of 22 Ne in the CO-core. 

The effective mass and composition of this pocket, as 
it would contribute to the merged white dwarf is not well 
constrained by the models we have available. For example, 
as the star leaves the AGB, the mass of this ls O pocket 
may be some 3.10 _4 Mq. Subsequent steady He-burning in 



the post-AGB phase (Fig. [4j| produces an O pocket in the 
pre-merger white dwarf of w 0.008 M0(FWHM) having a 
■ ls O abundance /3i a n -nn ) ~ 0.01. 



In the models of ISaio fc Jeffervl i|2002l ), this pocket is 
destroyed when the pocket is reheated by the post-merger 
accretion-driven helium flash; a new but smaller pocket is 
established at the outer edge of the newly established He- 
burning shell. However, if the 18 O pocket or material from 
deeper in the CO core is mixed during the merger, it may 
survive. Recalling that there is insufficient carbon in the 
post-AGB intershell to account for all the carbon observed 
in EHe and RCB stars, it is highly plausible that any car- 
bon mixed into the merger-product envelope from the CO 
core will be accompanied by 18 O. This will be diluted by ad- 
ditional 12 C and ls O if the mixing penetrates significantly 
beyond the 18 pocket. 

We have therefore incorporated the outer layers of the 
CO core more realistically. The distribution with mass of 
4 He, 12 C, 14 N, 16 0, 18 and 22 Ne is assumed to resemble 
that of the 0.60 Mq pre white dwarf shown in Fig. [4] where 
a core mass of 0.58 Mq is defined by the point where the car- 
bon and helium abundances are approximately equal. This 
composition distribution can be applied to COWD cores of 
different mass (mco:Co) by scaling its thickness inversely as 
(0.58/mco:co) 4 (this scaling also approximately reproduces 
the mHe:CO — n^cO:CO relation implied by Fig. [3}. 
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The abundances of N, 1 O and Ne may be scaled 
with metallicity Z, since 14 N derives mainly from the initial 
metallicity, and 1S and 22 Ne are mainly formed from a 
captures on 14 N. 

To explore the consequences of varying the 12 C: ls O ra- 
tio in the outer core, the factor fco can be used to force a 
rescaling of these two species. The model shown in Fig. [4] 
has fco = 0.8, but our AGB models suggest fco = 0.5, 
probably reflecting differences in the adopted 12 C(a,7) 16 
rate. 

A parameter a m i x represents the mass within the 
boundary layer which is mixed into the merged envelope. 
Formally, if j3(m r ) represents the distribution of mass frac- 
tion with mass inside the star (m r ), we compute 

Pi-.co-.co = / l3i:co:Co(m r )dm r / / dm r (33) 

for i corresponding to 12 C, 14 N, 16 0, 18 and 22 Ne. The 
mass limits are given by 

m D = m co:co + m s h, (34) 
m\ = m - a mix m sh , (35) 
m sh = (0.12/m CO :co) 4 . (36) 

The shell mass m s h characterises the scaled shell thickness 
corresponding to 0.12 M© in the 0.60 Mq pre- WD model. 

Thus Qmix = means that no carbon-enriched material 
from the boundary layer is mixed. a m i x = 1 implies that 
all material down to the point where the carbon and he- 
lium abundances are equal is mixed. a m ix = 2 means that 
the mixed layer reaches the region where 22 Ne and 16 are 
abundant. 

Other /3i : cO:CO are as given previously, except 

/JHe:CO:CO = 1 ~ £ift:CO:CO, (37) 

since this layer includes some helium from the base of the 
helium layer. The layer masses rrijk are adjusted to take the 
blurring of the carbon-helium layer boundary into account. 



4-2.4 Final abundances 

The ingredients of our model thus comprise five layers of ma- 
terial with masses rrijk, each having a representative com- 
position Pijk defined by current stellar evolution theory. As- 
suming our primary assertion that all of these layers are fully 
mixed during the merger, and that no further nucleosynthe- 
sis affects the apparent surface composition of the merged 
product, then the latter is simply represented by 

pi — 'EjkmjkPijk/'Ejkmjk (38) 

For comparison with observation, these abundances can 
be transformed to units more familiar in observational stud- 
ies. Recall that mass fraction P is defined in terms of number 
fraction n: 

Pi = nifii/Hinifii (39) 




N14 
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Figure 4. The distribution of various elements by mass fraction 
in the outer layers of a 0.6 Mq post-AGB star evolving towards 
the "knee" of the white dwarf sequence. 



whence 



ti = log m — log — + log Ei/Xim , (40) 

= log &± + C' (41) 

Mi 

C = log E,n lQMl - C (42) 

[X] = ei-e iQ , (43) 



where m are the relative abundances of species i by number 
(Erii = 1), and C is defined by Eq. [5] 

4.3 Nucleosynthesis during a hot merger 

In the cold merger model, the chemical composition of 
the HeWD is assumed to be unchanged during a merger 
with a COWD. However, SPH calculations indicate that 
some of this mate rial may briefly reach temperatures of 
6 x 10 8 K or more ^Guerrero et all 12004 lYoon et al.l 120071 ; 
lLoren-Aguilar et al.l I2009T ), and t hat some n u cleosy nthe- 
sis of a-rich material will occur. lYoon et all lj2007l l and 
lLoren-Aguilar et al.l (|2009l ) demonstrate how the disrupted 
material forms a relatively unprocessed disk containing 
slightly more than half of the HeWD, and a heavily pro- 
cessed corona containing the remainder. How the disk and 
corona are subsequently assimilated into the merged star, 
and what mixing processes occur, is not yet obvious. The 
simplest assumption is that turbulent mixing during the 
merger, and nuclear- and surface-driven convection follow- 
ing stable helium-shell ignition will completely mix both disk 
and corona with material from the intershell of the COWD. 

This process c a n be incorporated into our model. 
lLoren-Aguilar et al.l l|2009h (Tables 1 and 2) give sample 
masses and chemistries for the disk and corona in the cases 
of a 0.3 + O.5M and a 0.4 + 0.8M Q He+CO WD merger. 
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Using iiihciHc from Eq. [12] and interpolating we obtain 
mdisk:Hc and ro C orona:He for 0.3 < mHo:Hc < 0.4 (we do 
not extrapolate), and also chemistries for the same material. 
lLoren-Aguilar et al.l (|2009l ) compute models starting with a 
pure helium mixture for the HeWD and a carbon-oxygen 
mixture for the CO white dwarf, so we arbitrarily impose 
a Z-dependent lower limit on individual abundances in the 
disk and corona as given by the prescription for the cold 
merger. We note that contamination by a small amount of 
hydrogen and other metals will profoundly affect the pre- 
dicted hot merger nucleosynthesis, but we have no data on 
how this might affect current results. 



5 ELEMENTAL YIELDS 
5.1 Cold merger 

The predicted surface composition following the cold merger 
of a He+CO WD with some core mixing is illustrated in 
Figs. OH and El 

Fig. [5] shows the surface abundance of each species rel- 
ative to the solar surface abundance as a function of m ag b 
(the progenitor mass of the AGB star) for each of four ini- 
tial metallicities Z. This figure demonstrates that several 
species show high yields roughly independent of m ag b- 12 C 
and 18 O show uniformly high yields as defined by the core 
mixing parameters. The 16 yield is proportional to the 12 C 
yield; effectively determined by the le O/ 12 C ratio in the 



outer edge of the core. The N yield is proportional to Z 
as expected; it is determined by the initial CNO abundance. 
Note how the yields of 19 F, Na and Ne are strongly peaked 
in the interval 1.5 < m ag b/Mo < 3, whilst 31 P and Al 
are only significantly enhanced (> 0.5 dex) toward higher 
masses (m ag b > 2 M©). These yields reflect the AGB in- 
tershell nucleosynthesis. Meanwhile 12 C, 14 N and ls O are 
enhanced almost uniformly with m ag b and Z, effectively as 
the model was designed to deliver. 

Fig. [6] shows these same data rearranged as a function 
of initial metallicity [Fe] = log Z/Zq for four representative 
masses m ag b = 1, 1.9, 3 and 5Mq. They are plotted together 
with the observed abundance data in exactly the same way 
as in Fig. [1] Since we have only used two AGB models for 
the s-process elements, a single broken line represents the 
provisional prediction for m ag b ~ 3 Mq . Several correlations 
are noteworthy. 

5.2 Cold merger: individual elements 

Carbon: the model carbon enhancement was chosen to 
match the abundances measured in EHe stars by adjusting 
the degree of core mixing. The EHe carbon measurements 
are probably more representative than the RCB measure- 
ments because of the "carbon problem" referred to earlier. 
The EHe star data show a significant scatter. Allowing a m i x 
to vary by ±1 gives results still broadly consistent with the 
observed carbon (and oxygen) abundances. 
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Figure 6. As Fig. [T] The over-plotted lines represent the surface abundances predicted from a cold CO+He white dwarf merger. 
Predictions for four initial masses are shown: iriagb/Mg = 1 (dotted), 1.9 (dash-dot), 3 (dashed) and 5 (solid). Two or more lines are 
coincident in several panels, especially where only a single solid line appears. The very provisional result for s-process elements (see 
§4.2.2) is represented by a single dashed line. In this simulation, most of the the carbon and oxygen are dredged from the carbon/helium 
boundary layer at the top of CO core (a m i x = 3, fco = 0.8). 



Nitrogen: a generally excellent correlation with the observa- 
tions, supporting the basic assumption that the surfaces of 
EHe and RCB stars are primarily CNO-processed helium. 
Oxygen: a modest 18 O pocket at the helium-carbon bound- 
ary in the CO white dwarf progenitor may explain the most 
extreme oxygen abundances seen in metal-rich EHe and 



RCB stars, providing this layer is mixed during or immedi- 
ately after the merger. Note, however, that oxygen isotope 
ratios are known for only a few RCB stars and for no EHe 
stars, that ls O may be destroyed by high temperatures dur- 
ing the merger. An alternative is that the 16 0: 12 C ratio in 
the outer core is approximately unity or more. 
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Figure 7. As the top 4 panels of Fig. [T] repeated with different values for the core-shell boundary mixing parameter « m i x . The top row 
represents the inclusion of no helium-depleted material from the boundary layer. The second row represents mixing down to a point 
where carbon and helium abundances are equal. The bottom three rows represent increasingly deep mixing into material stratified as in 
Fig. PA 

^mix — 1 includes 50% of the pocket. The contribution of to the total oxygen abundance is shown in olive (grey). 



Neon: intermediate-mass models (1.9, 3 Mq) show significant 
enhancements of neon, formed primarily in the intershell of 
the AGB precursor. These match a few of the observed neon 
abundances, but the very high abundances measured in at 
least seven EHes are not yet explained by this model. 
Fluorine: as in the case of neon, a significant 19 F excess is 
generated in the intershell of the intermediate-mass models 



(1.9, 3 Mq), suggesting a likely source for the observed ex- 
cess. However, the predictions remain « 1 dex smaller than 
the measurements. However, ls O from the C/He boundary 
in the COWD provides an ample reservoir for the prompt 
production of 19 F if it is sufficiently heated, together with 
protons, during the merger; thus the model prediction rep- 
resents a strict lower limit. 
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Figure 8. As Fig. [6]cxccpt that the over-plotted lines represent the surface abundances predicted from a hot CO+Hc white dwarf merger 
and include no mixing with the CO core (a m i x = 0). Nucleosynthesis in the disr upted He- WD during the merger is interpolated between 
results for a 0.3 + 0.5Mq and a 0.4 + O.8M0 merger llLoren-Aguilar et alj|2009l) . The predictions for Ti and Cr are based on the model 
for a 0.3 + O.5M0 merger only. 



Sodium, Aluminium, and Magnesium: the recipe predicts 
some enhancement of these light elements over their initial 
values, particularly at \ow-Z for aluminium and magnesium. 
These predicted abundances are broadly consistent with the 
measurements for EHe and RCB stars in the case of sodium 
and aluminium. The under abundance of magnesium in both 
EHe and RCB stars at low-Z requires further investigation. 



Phosphorus: a primary motivation for this investigation, the 
recipe shows that 31 P generated in an AGB intershell can 
propagate to the surface of a subsequent white dwarf merger. 
The recipe only predicts significant overabundances at low- 
Z. Although this is consistent with one low-Z phospho- 
rus measurement and eight high-Z EHe measurements, the 
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recipe does it not explain a > 1 dex overabundances ob- 
served in eight other high-Z and one \ow-Z measurements. 
Silicon, Sulphur, Argon, Calcium, Titanium, Chromium, 
and Nickel: the cold merger recipe predicts negligible (sil- 
icon) or no enhancement beyond that expected from the 
enhancement of a elements in low-Z progenitors. The ob- 
servations are broadly consistent with the recipe predictions 
for all of these elements. A number of stars with up to 1 
dex enhancements of silicon, sulphur and/or argon demand 
further attention. In particular, the observed scatter of ±1 
dex in silicon is not obviously explained by this recipe. 
s-process: our provisional predictions for the s-process ele- 
ments zinc, yttrium, zirconium, barium, lanthanum, cerium 
and neodymium indicate a substantial excess is expected in 
all cases except zinc. The predictions are all in excess of the 
mean trend of the observed abundances (or their upper lim- 
its). For Y, Zr, and La, they are not in excess of the upper 
limit of the observed abundances. Given the small number 
of models available, and the probability that other factors 
will affect the observed distribution of s-process abundances, 
these results are inconclusive, but encouraging. 

It will be noted (e.g. from Fig. [6j that the predicted 
excesses of certain elements are strong functions of the pro- 
genitor mass (m ag b) and metallicity. Notable amongst these 
are neon, fluorine, sodium, and magnesium. An early objec- 
tive of this investigation was to determine whether the ob- 
served abundances placed any firm constraints on the pro- 
genitor mass. For example, significant excesses (which are 
observed) in neon, fluorine, and sodium are predicted for 
1.9 < m ag b/M0 < 3, whilst an excess of magnesium (which 
is not observed) is only predicted for m ag b / M© > 3. The ev- 
idence from phosphorus and aluminium remains ambiguous. 
The suggestion is therefore that the progenitor AGB stars 
had initial masses in the range 1.9 < ma g b/M© < 3. This 
suggestion assumes that the binary components evolved es- 
sentially as single stars up to the AGB, plausible if the 
first common-envelope phase required to produce the short- 
period double-white-dwarf binary occurred after the more 
massive star reached the AGB. 

5.3 Cold merger: core mixing 

Figure [7] illustrates the effect of a m ix on the abundances of 
the principal elements C, N, O and Ne, and also shows the 
contribution of 18 O to the total oxygen abundance. It must 
be noted that these results depend strongly on the compo- 
sition profile at the carbon-helium boundary, and in partic- 
ular on the ie O: 12 C ratio immediately below the boundary 
layer. The intent is to demonstrate the effect of increasing 
the depth of mixing on the final model abundances. 
Q?mix = 0: the first row of Fig. [JJ represents the case of no 
contribution from the boundary layer, oxygen completely 
fails to match the observed abundances of EHe and RCrB 
stars. There may be sufficient carbon in the intershell of 
intermediate mass (m ag b ~ 2 — 3 Mq ) to reproduce the lower 
envelope of the carbon and neon abundances, but not the 
full range. 

a m i x = 1: the second row of Fig.[JJ representing mixing down 
to the helium/carbon equilibrium point, has similar results 
for carbon, nitrogen and neon, but immediately shows much 
more oxygen at high Z. This oxygen is almost entirely 18 O 
(assuming that it survives the actual merger). The care- 



ful reader will note that the carbon abundances are slightly 
depressed compared with a m i x = 0; this is due to a mis- 
ma tch between th e intershell carbon abundances given in 
the lKarakas] (|2010h models and the abundances in the pre- 
WD model shown in Fig. 4. 

Q?mix = 2: slightly deeper mixing substantially improves the 
carbon result without much change to other elements. 
Q?mix >> 2: with very deep mixing, the models start to show 
too much carbon. Oxygen is increasingly composed of 16 O 
and becomes independent of Z, although at high Z, 18 O 
remains a significant constituent. The contribution of 22 Ne 
from the core has an almost negligible impact, even in high- 
Z models. 

As noted above, a shift in the 12 C/ 16 ratio immedi- 
ately below the boundary layer alters the chemical balance. 
By setting fco = 0.3, it was possible to obtain a good cor- 
respondence between model and observation for C, N, O, 
and the 18 0/ 16 ratio at the same time (not shown). How- 
ever, the value of fco required is unjust i fiably smaller than 
the value of 0.5 obtained in the lKarakas! (|2010l ') models. The 
prospects for neon seem less good; its abundance is primarily 
limited by the original CNO abundance. 

Further investigation of the chemical structure of white 
dwarf models derived from realistic AGB calculations will in- 
dicate whether the cold merger model stands up to scrutiny. 
Meanwhile, Fig. [JJ suggests that « m i x = 3 is satisfactory for 
the time being. Moreover, by indicating how much material 
from below the C/He boundary must have been mixed, the 
value of owx may also tell us something about the (cold) 
merger process itself. 



5.4 Hot merger 

Figure [5] shows the chemical yields predicted by this sim- 
plification of the hot merger model. In this case, there is no 
contribution to the carbon and oxygen from mixing with the 
He/CO boundary layer from the COWD. All of the excess 
carbon and oxygen comes from nucleosynthesis during the 
merger. 

The calculations bv lLoren-Aguilar et al.l (|2009l ) indicate 
significant production of iron, nickel and zinc, and very large 
quantities of argon, calcium, titanium and chromium pro- 
duced by a-capture reactions in the hot corona, especially 
towards the upper limit of the He- WD mass range (0.4 Mq). 
At this limit, the titanium and chromium yields are so high 
as to be off scale in Figure [8] thus we have restricted the nu- 
cleosynthesis of these two elements to the lower limit of their 
predicted range. The theoretical yields of iron, nickel and 
zinc are below the threshold defined by the initial metallic- 
ity and hence have no effect on the final abundances (Fig. [SJ 
compared with the model for the cold merger (Fig. |5J). 

The surface abundances following a hot merger as pre- 
dicted by the simple recipe are shown in Fig. [8] Recall that 
mixing at the carbon-helium boundary is switched off in 
this case. The predictions are complicated by having very 
few models amongst which to interpolate, making the Z- 
distribution of elements strongly affected by merger nucle- 
osynthesis indicative rather realistic (cf. argon and calcium) 
Carbon and Oxygen:, all surface carbon and oxygen is pro- 
duced by nucleosynthesis during the hot merger. The quan- 
tities are comparable with those from the cold merger, but 



18 C.S.Jeffery, A.I.Karakas & H.Saio 



in this case no parameters were tuned to achieve the correct 
outcome. 

Sulphur: additional sulphur produced in the hot merger pro- 
vides a better fit to the observations than in the cold merger. 
Argon: there is weak evidence of an argon excess in some 
EHes; this provides some support for its formation in a hot 
merger, possibly in higher-mass mergers (m ag b > 2Mq). 
Calcium, T itanium, and Chrom i um: the SPH merger cal- 
culations of lLoren- Aguilar et all (|2009l ), combined with the 
merger recipe described here, predict very high surface abun- 
dances of these three elements in nearly all cases. There is 
no observational evidence that any of these elements is sig- 
nificantly overabundant in any EHe or RCB star analysed 
to date. If EHes and RCBs are formed in a merger, there 
is no evidence that reactions leading to the production of 
calcium, titanium, or chromium operate during the merger 
process. This places strong constraints on temperatures and 
timescale of the hot merger. 

Fluorine: although not treated in the models of 
lLoren-Aguilar et al.l (|2009l ). any 18 O produced or mixed 
into the heated material will be at least partially burnt to 
make 19 F, providing the temperature does not significantly 
exceed 3.10 8 K, where the F will subsequently be destroyed 
by a captures. Hence the presence of any 18 O will almost 
inevitably result in an excess of 19 F. 

Lithium: Li should be d estroyed duuing a hot merger. 
lLoren-Aguilar et alJ l|2009l ) do not include Li, so predictions 
for Li are not shown in Fig. [8] 

5.5 Uncertainties 

The obvious question is what confidence can be placed on 
using a simple recipe to predict the outcome of a compli- 
cated process? For any element which is produced outside 
the standard hydrogen and helium-burning reactions (e.g. s- 
process elements in the AGB precursor), the major impact 
on the final surface yield is the degree of dilution by the 
accreted helium white dwarf. In our recipe this is primarily 
constrained by the minimum white dwarf mass for a conser- 
vative merger. Were mue to exceed this value, the predicted 
excesses would be reduced. This does not help to explain 
neon, fluorine or phosphorus in EHes. Were the merger to 
be non-conservative (efficiency a < 1), yields might increase 
by factors 1/a. SPH calculations support a = 1. Efficiencies 
a < 1 would be surprising. Explaining neon, fluorine and 
phosphorus would require a < 0.1. 

The neglect of specific mixing processes or individual 
reactions in the stellar evolution models will impact on the 
recip e predict i ons. F or example, the intershsell compositions 
from lKarakas! (|2010h were computed without the addition of 
a 13 C pocket. This pocket is thought to form by the mixing 
of protons from the H-rich envelope into the intershell dur- 
ing the deepest extent of a third dredge-up episode. These 
protons are quickly captured by the abundant 12 C resulting 
in the formation of a 13 C pocket. Here, neutrons are lib- 
erated by the reaction 13 C(a,n) 16 during the interpulse. 
Note that in most calculations the 13 C pocket is added arti- 
ficially or i nduced throu gh the inclusion of convective over- 
shoot (e.a. \Herwicl\200(\) . Whilst 13 C pockets facilitate the 
formation of s-process elements, they are also important for 
enhancing the abund ance of some lighter elements inclu ding 
19 F, 22 Ne, and 23 Na (|Lugaro et al. I l2004l ; lKarakasll2010l ). 



Legitimate questions concern the core and shell masses 
for the AGB stars, the yields obtained in the AGB inter- 
shell nucleosynthesis and in the hot merger nucleosynthesis. 
For example, how sensitive are the masses adopted here to 
the microphysics, e.g. convection, rotation, mass-loss and/or 
mass transfer? Since the model abundances are computed 
for CO+He WDs with q = q a it, would additional dilution 
produced q > g cr i t seriously compromise the results? If mass 
transfer occurs before AGB evolution is complete, or even 
before the star reaches the AGB, are the intershell-yield ver- 
sus initial-mass relations still acceptable? What would the 
hot-merger nuclear-hydro calculations look like with more 
realistic networks and initial composition? We have dis- 
cussed the question of the chemical structure of COWDs, 
particularly at the C/He boundary layer. Additional calcula- 
tions for hot mergers with more extensive reaction networks 
and robust starting mixtures are also urgently needed. 



6 CONCLUSION 

The object of this paper was to clarify the surface abun- 
dances which might arise under conservative assumptions 
for the merger of a helium white dwarf with a carbon-oxygen 
white dwarf, and to compare these with observed abundance 
anomalies for extreme helium and RCoronae Borealis stars. 
We have collated the observational data describing the sur- 
face abundances for the latter, and presented it, element 
by element, in a way that demonstrates any current ex- 
cess over the progenitor composition. We have discussed the 
background to the physics of binary white dwarf mergers 
and their association with the formation of extreme helium 
stars and RCoronae Borealis stars. We have developed a 
more elaborate version of the simple recipe used to estimate 
surface abundances in previous discussio ns of this question 
l|Saio fc Jeffery||2002l ; rPandev et al.ll2006l '). We have incorpo- 
rated state-of-the-art calculations of the masses a nd light- 
eleme nt composition of AGB intershell regions (|Karakasl 
l2010l) . We have made allowance for the existence of an 18 O 
pocket at the outer edge of the CO core. We have considered 
the difference between a cold and a hot merger, i.e. whether 
additional nucleosynthesis occurs during the destruction of 
the helium white dwarf. 

Both models successfully match, or can be made to 
match, the observed surface abundances of carbon, nitrogen, 
and oxygen. The excess nitrogen comes primarily from the 
helium white dwarf as the residue of CNO-processed carbon 
and oxygen. In the case of the cold merger model, the excess 
carbon and oxygen comes from the carbon-helium boundary 
of the carbon-oxygen white dwarf. A substantial fraction of 
oxygen probably takes the form of 18 O from a pocket just 
beneath this boundary, but 16 O dredged from deeper layers 
will also be present. In the case of the hot merger model, 
the excess carbon and oxygen can be produced during the 
merger. It is possible that the observed excess comes from 
both sources. 

Both models predict up to 1 dex enhancements of 19 F 
and 31 P, but not enough to match the observed overabun- 
dances of fluorine and phosphorus. These elements come 
from the AGB intershell. An examination of mixing at the 
carbon/helium boundary of the CO WD suggests that ob- 
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served neon may come from the outer part of the car- 
bon/oxygen core. 

Both models predict modest overabundances of sodium, 
aluminium and magnesium, particularly at low metallicity. 
These broadly match the observed abundance distributions 
in sodium and aluminium. Magnesium is not observed in 
excess at low-Z. 

The hot merger model currently predicts overabun- 
dances of calcium, titanium and chromium which are not 
observed. The model overabundance of argon is roughly con- 
sistent with available measurements, suggesting that argon 
might be produced during hot merger nucleosynthesis. We 
do not yet have AGB intershell yields for argon, which would 
affect the cold merger predictions. 

We still require state-of-the-art data for s-process yields 
in AGB intershell. 

Overall, the majority of species observed to be over- 
abundant in EHe and RCB stars are found to be enhanced 
in reasonable quantity in one or other of the white dwarfs 
prior to merger. In other words, additional nucleosynthesis 
during a merger solves very few problems, although that 
does not mean that it does not happen. 
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Table 1. Observed abundances for EHe and RCrB stars (omitting 13 C/ 12 C, 18 0/ 16 0, Sc, Mn, Co, Cu, and Sr). 



